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Abstract 

We analyze the effect of R-parity violation in the minimal supersymmetric 
standard model on the CP asymmetries in b — > S7 decay. The direct and 
mixing-induced CP asymmetries arising from the lepton number violating 
couplings are strongly constrained by the current experimental limits on the 
corresponding couplings. Allowing a heavy neutrino (m UT ~ 10 keV) and a 
moderate mass splitting of sfermions, the direct CP asymmetry around 15 
% and the nearly maximal mixing-induced CP asymmetry (~ 100%) can be 
realized, depending on the R-parity conserving contributions to the radiative b 
decay. With the baryon number violating couplings, only the mixing- induced 
CP asymmetry arises and it can be maximal provided there is a the similar 
sfermion mass splitting. 

PACS number(s):11.30.Er,12.60.Jv,13.20.He,13.40.Hq 
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I. INTRODUCTION 



CP violation in radiative B decays may provide a promising tool for probing new physics 
beyond the Standard Model (SM). The direct CP asymmetry in the radiative b — > 37 decay 
defined by M 



T{B -> X sl ) - F(B - Xyy) 



T(B^X sl ) + T(B^Xrf) 



(1) 



£ 7 >(l-<5)E? ax 



is below 1 % in the SM 0. Therefore, the observation of a sizable CP asymmetry would 
be a clean signal of new physics and may further discriminate various extensions of the SM. 
Recent model-independent analyses of CP violating effects in the inclusive decays B — > X s ^f 
in terms of the Wilson coefficients of the dipole moments operators have shown that models 
with enhanced chromo-magnetic dipole moments can naturally provide a large asymmetry 
U and it can reach up to 30 % accommodating the observed branching ratios of b — > 37 
and b — > sg decays The specific predictions for the CP asymmetry are of course model- 
dependent. The left-right model can yield A C p at the level of 1 % |4[], and a two Higgs 
doublet model up to 10 % M. In supersymmetric extensions of the Standard Model, there 
exists additional CP phases which would give rise to large CP violating effects. In the 
context of minimal supergravity models ||, the direct CP asymmetry turns out to be less 
than 2 % due to strong constraints on supersymmetric CP violating phases from the neutron 
and electron electric dipole moments [[/J. The asymmetry can be enlarged in non-minimal 
models up to 7 % with more freedom in CP phases M, or up to 15 % with generic sfermion 
mixing ||||. 

Another important CP violating observable is the mixing- induced CP asymmetry in 
exclusive radiative B decays |ITJ which can occur for radiative B q — > M g 7 decays, where 
M q= ds is any hadronic self-conjugate state with CP eigenvalue £ = ±1. The CP asymmetry 



in the time-dependent decay rates Y(i) for B q 
T(t) - f(t) 



M q >y and T(t) for B q -> M g 7 is then 



A(t) 



r(t) + f(t) 



£A M sin( 



>M 



l - <Pr) sin(Amt) , A 



M 









l 2 + 
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(2) 



where 0m and Am are the phase and the mass difference of B q -B q mixing, respectively, 
and Cjljr are the effective coefficients of the left-handed and right-handed dipole moment 
operators for the b — > 57 decays, and 4>l,r are their phases, respectively. In the standard 
model, such asymmetries are expected to be small as one has Cjr/Cjl ~ rriq/rrib and thus 
Am of order 1 % and 10 % for the b — > and b — > decays, respectively. Unlike the direct 
CP asymmetry, the left-right symmetric model may allow the mixing-induced asymmetry 
up to 50 % [llj. The supersymmetric model with generic sfermion mixing can yield even 
larger asymmetry up to 90 % 0. 

In this paper, we will analyze the effects of R-parity violation in the Minimal Supersym- 
metric Standard Model (MSSM) on the CP asymmetries in radiative B decays. R-parity 
violation in the MSSM introduces a large number of trilinear couplings which violate lepton 
and baryon number. These additional couplings can surely be sources of flavor and CP 
violation, which might lead to a huge effect on the CP-odd observables in the B decays. 
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CP violating effects of R-parity violation in the hadronic B decays have been considered 
previously in Refs. [[l I -15|, showing that significant modifications to the SM predictions can 
follow from R-parity violation. This work is devoted to the analysis of the CP asymmetries 
in the radiative B decay. As we will see, contrary to the cases with hadronic B decays, 
various experimental constraints on the R-parity violating couplings coming particularly 
from rare B decays strongly limit the amount of the direct CP asymmetry in the b — *■ sj 
decay whereas a large mixing-induced CP asymmetry can be induced from, in particular, 
baryon number violation. 

This paper is organized as follows. In section II, we analyze the general effective Hamilto- 
nian describing the b —>■ s'j decay including the new operators induced by R-parity violation. 
We derive the anomalous dimensions and evolution matrix with the enlarged operator set. In 
section III, we discuss the CP asymmetries arising from lepton number violating couplings 
with which we can have both the direct and mixing-induced CP asymmetry. We deal with 
both cases in separated subsections. In section IV, we discuss the CP asymmetries arising 
from baryon number violating couplings, in which case only mixing- induced CP asymmetry 
can be obtained. We conclude in section V. 



II. R PARITY VIOLATION AND EFFECTIVE HAMILTONIAN 

Let us begin our discussion with defining our choice of the basis to describe the R-parity 
violating couplings. For comparison with experiments, it is convenient to work with R-parity 
violating couplings defined in the quark and lepton mass eigenbasis. In this prescription, we 
leave the neutrinos in the charged lepton mass eigenbasis since neutrinos can be taken to 
be massless for our purpose. The full superpotential of the MSSM fields including generic 
R-parity violating couplings is then given by 

W = /i#i# 2 + KHxLiEt + hUHlDiDt - H^V^D*) + h^U^ - H+V^U?) 
+ - L]E % )El + KjkiLiDjDt - E^Dl) + l^U'Dpi , (3) 

where Li = Ei) and Qi = (Ui, Di) are the lepton and quark SU{2) doublets, and Ef, U°, 
D\ are the SU(2) singlet anti-lepton and anti-quark superfields. Here is the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix of quark fields. In order to ensure the longevity of 
proton, the products A' A" have to be highly suppressed fl"6|l . For this reason, one usually 



assumes lepton or baryon number conservation to discard the couplings A/A' or A", respec- 
tively. In this paper, we discuss both cases separately: one with R-parity and lepton number 
violation with nonvanishing A', and the other with R-parity and baryon number violation 
with A". The couplings A will be irrelevant for our discussion. 

The presence of R-parity violating couplings in Eq. (§) give rise to new contributions to 
the radiative decay b — > 57 through one-loop diagrams exchanging sleptons or squarks []I7| as 
depicted in FIG. 1. Note that R-parity violation may induce equally sizable dipole moments 
of the left-handed and right-handed type, respectively labeled by L and R as follows: 

7L cxiga^F^, 7R K^%a^b a L F^. (4) 
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In the SM and many extensions of it, the coefficients of the second one is usually suppressed 
by the factor m s /m b . In our case with lepton or baryon number violation in the MSSM, the 
operator Ojl and O^r are generated for nonvanishing combinations of couplings; 

^n3j^n2j aI1 d ^njV^njzl ^nl2^n23 J (5) 

respectively. With the couplings in Eq. (^j), there arise also other four-quark operators 
which should be taken into account in the complete effective Hamiltonian describing the 
b — > sj decay. The whole set of the effective operators arising from R-parity violation can 
be described by a simple generalization of the SM operator space by separating the standard 
operators 03,4,5,6 for each quark flavor q — u,c, d, s, b. That is, we introduce the additional 
operators, 

Ok = ^bliri.ql , 0\ L = ~%fV> qj^qi , 

<nL = ~3ZfVL<&idk> °Il = sh^llbA, (6) 

where a, f3 are the color indices. The dipole moment operators for the b — > 57 and b —>■ sg 
are defined by 

We also have the right-handed counterpart of the operators which can be obtained by the 
exchange L <-> R in Eqs. (^) and (|7|). 

The effective Hamiltonian at scale fi < 0(mw) relevant for the b —>■ 57 decay is now 
described in terms of the enlarged set of operators as follows: 



E CM om + EE c]M o]M + {l^r) 

i=l q j=3 



(8) 



where A* = V t * s Vtb, Gf is the Fermi constant and C"s are the Wilson coefficients which will 
be determined later. The operators On with i = 1, • • • , 8 are those considered in the SM 
1~8| . Note that On,2L = 0^ L 4L and On = J2 q 0\ L for % = 3, • • • 6. 

Given the Wilson coefficients Cj, C| including the contributions from R-parity violation 
at the weak scale fi = mw, we need to calculate those at the scale \i ~ m b through the 
renormalization group (RG) evolution. Since the QCD running do not mix the left-handed 
and right-handed set of operators and its effect is identical, it is enough to calculate the RG 
equation at the one sector. At the leading order, it is rather straightforward to calculate 
the anomalous dimension matrix in the extended operator basis following the standard 
calculation [0. At this point, let us recall that it is convenient to use the so-called "effective 



coefficients" |19|] which are free from the regularization-scheme dependency in the mixing 

•^7,8 



between the sets O^ with i — 1, • • • , 6 and O78 resulting from two-loop diagrams. In terms 



of the effective coefficients, 
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CV e// (/i) = CV(^) + 5>C7(/i) 

C e 8 ff (fi) = Cs^)+J2^Ci^), (9) 

with the index / running for 26 indices labeled by % and jq for the operators Oi (i ^ 7, 8) 
and Oj, the effective anomalous dimension matrix is given by 

r 7J7 + Ex VkIjk - VJ177 ~ z,n% 7 , for J = 7, J 7^ 7, 8 
7j/ J = S 7J8 + Ex ^x7jx - ZJ"T88 , for J = 8, J ^ 7, 8 (10) 
I 7 j/ , otherwise. 

The corresponding RG equations for the effective coefficients are then 

d -cfM = ^zil/cf^) • (ii) 



<iln/i 47i 



In the naive dimensional regularization scheme which we follow p9| , the nonvanishing coef- 
ficients yi, zj are 

1 

2/5 = 2/56 = — g , 2/6 = J/66 = -1 

z 5 = z 5b = -1 . (12) 
Now, keeping track of the flavor structure of the standard calculation of the matrix 



T8fl , one can find the 28 x 28 anomalous dimension matrix in a straightforward way. The 



results are presented in the Appendix. There, we also calculate the evolution matrix from 
which the Wilson coefficients at the scale \x = nib can be obtained in terms of the coefficients 
determined at the weak scale. The contributions to these coefficients from R-parity violation 
will be discussed in detail in the following sections. 



III. CP ASYMMETRIES WITH LEPTON NUMBER VIOLATION 

A. A^ 3 jA^ 2 j : Direct CP asymmetry 

With the nonvanishing combination of the R-parity violating couplings A^ 3j A^ 2 ,-, the 

effective Hamiltonian in Eq. (§) includes only the operators of the left-handed type: O d ^ L 
induced by the tree-level diagram exchanging sneutrino v n . In addition to these, there arises 
also additional effective semileptonic operators, 

°9L = Sll^L VnL^^VnL , (13) 

through the tree diagram exchanging the right-handed down squark dj. The corresponding 
Wilson coefficients are given by 

i 3 V V* 

™^ )= T?kz£^- (14) 



5 



Note that C 5 3 L (mw) = and it remains vanishing at low energy scale under one-loop RG evo- 
lution as can be seen from Eq. (A9) in the Appendix. The R-parity violating contributions 
to the coefficients of the operators 7 l and Osl come from one- loop diagrams exchanging 





sneutrinos and right-handed down squarks, and have been computed in Ref. 
be conveniently re-written in terms of the coefficients C 6 J L and Cg£ as follows: 



Cjl {m w ) 



QdCs£ {m w ) 



Qd 
6 



J2C» 9 l(m w ) + 2j2C d 6 i(m w ) 



They can 



(15) 



TV 

where Qd = — 1/3. In deriving C 7 £ SL , we neglected the down- type squark mixing and the 
masses of down-type quarks compared to the mass of the sneutrino. As shown in Eq. (|12|). 
the effective coefficient C 7I f at /i = mw gets a nontrivial contribution from C 6 £ with j = 3 
and thus we have 



C 7 {J{m w ) = C 7L {m w ) - Cl L {m w ) 
C8L f ( m w) = C 8L {m w ) . 



(16) 



Making use of the relation (|15|) and the formula (|A10 ) in the Appendix, one can obtain the 
Wilson coefficients at the scale /i = as follows; 

C${m b ) = 0.67 CffVty) + 0.092 C^ u (m w ) - 0.17 C 2L {m w ) 

- O.U[C* L (m w ) + C s eL (m w )} - 0.80C b 6L (m w ) - 0.022 C&(m w ) , 

C e / L f {m b ) = 0.70 CH M {m w ) - 0.080 C 2L {m w ) 

+ 0A2[Ci L (m w ) + C s eL (m w )] + 0.WC b eL (m w ) + 0.12 C&(m w ) , (17) 

where C 7 f^ L (m\y) contain the contributions from the R-parity conserving MSSM sector as 
well as from the SM one. We assume that C2L{ m w) comes solely from the SM. To quantify 
C 7 fJf L (mw) for our purpose, we introduce the parameters 777,8 which are defined by 



C^ u (m w )^r, 7 C^(m w ) 



(18) 



The parameters 777^ are complex in general. 

Given the Wilson coefficients in Eq. ([17]) ? we are ready to analyze the CP violating 
effects from R-parity violation in the radiative B decay. Referring to the work by Kagan 
and Neubert Q for details, the direct CP asymmetry Aqp and the branching ratios for the 
decays b — > 57 and b — > sg are given by Q 



A 



CP 



{1.23 %[C 2L C* 7L \ - 9.52 %[C 8L C* 7L \ + 0.10 3[C 2L C 8 * L ]} (%) 



7L 



B(B -> X sl ) a 2.57 x 10- 3 iW<5) 
B(S^X sa )«0.96|C 8i | 2 B semi , 



(19) 



Here, we neglect the R-parity violating contributions to Aqp through terms such as ^s[C^ p C 7L \- 
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where the coefficients C's without arguments are understood to be the effective ones evalu- 
ated at the scale m h . The quantity K^ LO (6) — \Cj\ 2 + 0(a s , l/mf) contains the corrections 
to the leading-order result and the specific forms of the corrections can be found in Ref . [|TJ . 
We will take 5 = 0.3 and B semi = 10.5 % for the branching ratio of the semileptonic de- 
cay B — > X c e V. In this work, we take the following values for the SM predictions for the 
C||f 7i 8L at the mw scale; 

C 2 s f(m w )~1.0, C™(m w )^ -0.20, Cff (m w ) « -0.10 . (20) 

The above choice of parameters yields the values at the scale \i = = 4.8 GeV; 

Clf (m 6 ) « 1.11, Cff (m 6 ) « -0.32 , Cff (m 6 ) « -0.15 . (21) 



Considering the above values of the SM coefficients, Eq. (|17D suggests that a significant 
contribution from R-parity violation to the b — > S7 decay can arise for \Cg r l(mw)\ ~ 10, 
\C 6 ^(mw)\ ~ 2, or \Cl L {m w )\ ~ 0.3. Furthermore, as we will show, if a sizable \C\ L \ ~ 0.3 
is allowed, one can get the direct CP asymmetry of the order \Acp\ ~ 10 % satisfying the 
observed branching ratios of B — ► X s 7 and B — > X sg decay. 

In order to figure out how large CP asymmetry can come from R-parity violation, let 
us consider the experimental bounds on the new Wilson coefficients appeared in Eq. (|T7j). 
First of all, those coefficients will be constrained by the experimental data for the branching 
ratios in Eq. (|i"9"[ ). In this work, we use the CLEO data: 

B{B -> X sl ) = (3.15 ± 0.35 stat ± 0.32 stat ± 0.26 modcl ) x 10" 4 , [20] 
B(B -> X sg ) < 6.8 % (90%C.L.) . 



More important constraints on the coefficients Cg£ and Cg£ come from experimental data 
on the various B meson decays. Let us discuss the relevant bound for each coefficient. 

First, the coefficient Cg L is constrained by the B decay mode B° — > K°tt°, whose branch- 



ing ratio is observed to be [22 



B(B° -> K°n°) < 4.1 x 10" 5 



Following the similar method used in Ref. [23], we estimate the matrix element of the 
operator 0\ L as 

< K«AOLW >« i 0( m l {m \- ml) j KF^iml) . (22) 

2(m s + m d ){m b - m d ) 

Taking f K = 0.16 GeV, F^(m 2 K ) = 0.33, we obtain 



6L 



<0.17. (23) 



Second, the coefficient C$ L is constrained by considering the B decay mode B° —> (pK° 
The experimental limits on the branching ratio of this decay mode is |2l 

B(W -f 0K°) < 3.1 x 10~ 5 . 



7 



From our estimation of the matrix element of the operator Oq L ; 



< <j>K°\0' 6L \B° >^ j^Um^e- (P B + P K ) F( 



B - K (mi: 



AN 

where e is a polarization vector of </> and we will take N = 3. With 



B^K U 



[mV) = 0.38, we obtain following limit 



IQlI<0.23 



(24) 
0.23 GeV, 

(25) 



The above bounds in Eqs. (|23|) and fl25[ ) tell us that the coefficients C^l cannot play any 
important role in the radiative B decays as can be seen from Eq. (|H]) 0. 

Now let us consider the constraint on C\ L . The most useful bound comes indirectly from 
the consideration of the decay mode B° — > X s u n V^ whose branching ratio can be calculated 
as [p5fl 



B(B° -> X 8 v n V$ 



A* 



cb 



I riv n 1 2 



fps( m c/ m b) 



(26) 



According to the analyses in Ref . |25| , one obtains indirect experimental information on the 
above branching ratio: 

B(5° -> X s v n v^) < 3.9 x 10" 4 . 
Taking this value and |A t /Vd,| = 0.976, fps( r n>c/ rn b) = 0.5, we put the bound, 



I n Vn I 

\ U 9L\ 



< 0.044 . 



(27) 



Thus the contribution of the coefficient Cg£ to Eq. ( |TTD can also be neglected. Now, under the 
condition that the bound (^7]) is applied to each component of the coefficient [see Eq. (0)], 
we get 



\C b 6L \ < 0.044 



m 



mt 



(28) 



for each n = 1,2,3. Here we remark that under the condition fl2"T|), the R-parity violating 
contribution to the semileptonic decay B — > XJ n V^ through the effective operator 



n 



eff 



A n3j A 

2m 



n3j / ^n2j ur _ 

2 c L ^b L e nL 7 At i/ nL , 



(29) 



can be made small enough to satisfy the direct experimental bounds [26 



2 The bounds in Eqs. (|23| ) and (|25|) are at the scale m?,. Practically, the Wilson coefficients induced 
by the R-parity violating couplings at the scale m(, are nearly the same as those at the scale mw 
[see Eq. (A9)]. 
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lA^y < 1.1 x 10- 



"3 



100 GeV 



l«l <44xir3 (lM%v) 2 - (30) 

Finally, we have to consider the neutrino mass coming from our choice of nonvanishing 
A^ 33 A^ 23 . With nonzero value of A^ 33 , the neutrino v n may get an undesirably large mass 
from one- loop diagrams with the exchange of bottom quark and squark |27| . The one- loop 
contribution to the neutrino mass is given by 

m „ 3 \%Aml 

where A denotes left-right sbottom mixing parameter and fh is an average sbottom mass. 
Taking A'^ 33 /m 2 < 10 _7 GeV -2 from the consideration of Eq. (0), we obtain 

m, n <10keV (31) 

with A = 100 GeV. Thus it is well below the direct experimental limit for the muon (tau) 



neutrino which is 0.17 (18) MeV p2[ . Yet another indirect limit for the muon or tau neutrino 
mass comes from the observation of atmospheric muon neutrinos. The atmospheric neutrino 
data from the Super-Kamiokande are known to be nicely explained by the neutrino oscillation 
between largely mixed muon and tau neutrinos |2"%fl . A natural consequence of this would 
be that the muon and tau neutrinos are very light m v „ T < 1 eV. If this is the case, there is 
no room at all for large CP asymmetry from R-parity violation. However, having a large 
R-parity violation and thus a heavier muon or tau neutrino is not excluded completely as 
there exist some other viable options for the explanation of the Super-Kamiokande data, 
such as a neutrino decay |29 . 



From the above consideration, the only possibility for a significant enhancement of R- 
parity violating contribution to the b — > S7 decay is to have a large coefficient C\ L with a 
sfermion mass hierarchy mj c > rn^. For example, we need mj c ~ 5mj> n to get C\ L ~ 1. 
Taking into account all the above experimental limits on the Wilson coefficients, let us now 
analyze how large CP asymmetry Acp can be obtained. In FIG. 2, we show Acp as a 
function of the branching ratio of the decay B — > Xg'-f varying \C\ L \ and Arg(C| L ) from 
to 2n. The other R-parity violating couplings are neglected. We take r] 7 = r] S = r] a.s 
a real number in FIG. 2 even though 77 is a complex number generally. In this figure we 
consider additional experimental constraints coming from B(B — > X s ^y), B(B — > X sg ). If 
m| c are larger than ?n~ n , then sizable \C\ L \ is allowed evading the bounds Eq. (|28|) . We find 
the CP asymmetry can reach 13 % for vng c jmy n pa 3.4 with vanishing R-parity conserving 
supersymmetric contributions (rj = 1) as shown in the left-upper frame of FIG. 2. In the case 
where the R-parity conserving supersymmetric contributions take the same sign as the SM 
values of CV^sl^w), this CP asymmetry can be larger as shown in FIG. 2 with 7] = 2. On 
the other hand, the CP asymmetry decreases when the R-parity conserving supersymmetric 
contributions take the opposite sign to the SM values of C 7 L i8 L( m w) as seen from FIG. 2 
with 77 = 0, — 1. 

Before concluding this subsection, it is worthwhile to notice that the CP asymmetry in 
the hadronic B decays such as B d — > ttK s and B d — > (f>K s can be significantly affected by 
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the R-parity violating couplings Cfl even if the effects of the couplings on the direct CP 
asymmetry in the radiative B decay are negeligible |]T3|-|H)[] . 



B. A^j 2 A^j 3 : Mixing-induced CP asymmetry 

Contrary to the previous case, the combination of R-parity violating couplings AL- 2 A£,- 3 
leads only to the right-handed set of operators: 7R:8R and 0\ R&R in the effective Hamil- 
tonian (H). In addition, the coefficient of the operator 0\ R is also generated. Even though 
it does not contribute to the b decay, it's coefficient will be considered since it enters into 
Cjr^r- As in the previous subsection, we have also the semileptonic four-Fermi operators 
as follows: 



0\lR = ^b a R e^l lii e nL . 
The nonzero Wilson coefficients at my/ induced from our R-parity violation are 

-1 



(32) 



Cm{m w ) 



3 y y* 

\ ^ A nj2 A nj3 

4 v / 2G F A i tx m f 



CZ n R {mw) = C% R (m w ) 



1 



A\/2Gf At j=1 



3 y y* 

y-v A nj2 A nj3 



m, 



(33) 



where qj can be either Uj or dj, m~ ln and mQ, are the masses of the doublet slepton and 
squark, respectively. In deriving these coefficients, we do not consider the effect of the CKM 
mixing. Note that there can be also other four-Fermi operators involving two different 
flavors of quarks or leptons through the CKM mixing, which we neglect in our discussion as 
they give subleading contributions. Now, disregarding the contributions from the R-parity 



conserving supersymmetric sector, the Wilson coefficients for O-jr^r at the scale mw JT7 
can be expressed in terms of the coefficients in Eq. (|33| ) as follows 



C 7R (m w ) = Cfy (m w ) = % 

o 



1 



E C m(mw)+2j2C q 6R (m w ) 



■ £ C1or(™w) + 7 £ C q 6R (m w ) P 



(34) 



Cm{m w ) = C 8 £ (m w ) = r 



+ E CZ5R(m w ) + 2 £ C q 6R (m w ) P a ( Xj ] 



(35) 



where 



Pry (Xj j 



36 



[Q u Fi{x j ) + F 2 {x j )} 
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with Xj = ml /m~ . The functions F 1 2 are defined as 

1 (2 + 3x — 6x 2 + x 3 + Qx In x) 



F 2 (x) 



12 (1-x) 4 
1 (1 — Qx + 3x 2 + 2a; 3 — 6x 2 In x) 
12 fl - xY 



Note that these functions are defined as -P 7 (0) = P g (0) = 1 and their values for some selected 
x are listed in TABLE II 

In the case under consideration, there is no phase in C 7R C 8R contributing the direct CP 
asymmetry in right-handed sector analogue of Eq. (|T^) since the same combination of the 
couplings generates both C 7Rt8R . However, there can arise a sizable mixing- induced CP 
asymmetry as defined in Eq. @. That is, we may have \Am\ ~ 1 with \C 7 l\ ~ |CVr|- Here, 
the Wilson coefficients are effective ones evaluated at the scale m^. Combining our results 
in the Appendix (|A11|) and Eq. ([34]), we find 

C e 7 f R f (m b ) = 0.40 [Cf R {m w ) + Cl 2 R (m w )} + [-0.80 + 0.26 P 7 (x t ) + 0.031 P g {x t )\ Cf R {m w ) 
-0.022 C% R (m w ) + 0.31 C^ R {m w ) . (36) 

From this equation, it appears that we can easily obtain C 7R ~ \C 7 l\ = \C 7 f M \ ~ 0.32 \rj 7 \. 
To clarify this, we now consider the experimental constraints on the coefficients appearing 
in Eq. (j36|). Let us first note that the arguments in the previous subsection are applied here 
to get the bounds; 

|C&| < 0.17, |C£| < 0.23, \C&\ < 0.044 (37) 

as in Eqs. (|23|), (^) and (0), respectively. Thus, the contributions of these coefficients in 
Eq. (RBT) are not significant. Another important constraint on the relevant R-parity violating 



couplings comes from the decay B — > XJ+l n |Mj induced by the second effective operator 
in Eq. (|32]) . This consideration leads to 



where B(6 — > X s l^l n ) cxpt denotes the experimental upper limit on the branching ratio of the 
B — > XJ^l~ decay modes given by 



B(B -> X s e + e~) < 5.7 x 10~ 5 , 
B(B -> X s fi + fi-) < 5.8 x 10~ 5 . 

Even if there are no data for b — > X s t + t~ , we expect it's branching ratio is most probably 
less than that of B — > X s e + e~ . Then, similarly to Eq. (^), the bounds on C\ R can be 
obtained indirectly from Eq. ( j3~8l) as 



IC&I < 0-017 | \ r B ^'^ m I I ^ I . (39) 
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for each j and n. Therefore, we may obtain a sizable C^ R (nib) from the contribution 
of Csft(mw) if there is again a hierarchy between sfermion masses. For example, taking 
/thj ~ 5, |Cfjj| = 0.33 is allowed within the present experimental bound. Taking 
= 1 which gives P 7 (l) = 5/14 and P g (l) = 1/4, one could obtain \CfJ(m h )\ « 0.23. If 
0.23, \Am\ ~ 1 is possible accommodating the measured B(b — ► 57). Finally, 
let us note that a large mixing- induced CP asymmetry requires the (tau) neutrino to be 
heavy as discussed in the previous subsection. 



m 



\C e 7 l f (m b ) 



Again we note that with the coupling A^ 22 A^ 32 which does not affect the radiative B 
decays, one can have important effects on the CP asymmetries in the hadronic B decays 

♦ 7T ± K° [T5f. 



such as Bd — > 4>Ks pT5[ , P^] and B 



IV. CP ASYMMETRIES WITH BARYON NUMBER VIOLATION 

A„ 12 Ani3 : Mixing-induced CP asymmetry 

Our final case is to have baryon number violation while lepton number is conserved. 
Then, the new operator set for the b — > s transition contains again only right-handed ones 
with nonvanishing product of couplings A" 12 A^ 3 : 0^ R d iR . The Wilson coefficients of these 
operators calculated at the weak scale m w are 



C%k(m w ) 



C d R {m w ) 



\ll XII* 

A nl2 A nl3 



4V2G F \ t m% 



3 \ 11 \ a* 

A nl2 A nl3 



4y/2G F \ t 



E 

n=l 



AR ■ 



(40) 



Notice that the simultaneous presence of nonvanishing coefficients C^R,m is due to color 
antisymmetry in the superpotential term, U C D C D C . Following the similar steps as before, 
we get the relation 



C m {m w ) = Cj| (m w ) = - 



Cm{m w ) = C SR } (m w ) 



ACt R (m w ) -5^2C^(m w )P^(x n ) 

n 

C d m {m w ) + Y,Cr R {m w )P' g {x n 



(41) 



where 



PLiXr 



36 



[QuF^Xn) - Q d F 2 (x n )} , P'(x n ) = 12[Fx(x n ) - F 2 (x n )} 



with x n = m\Jmg c and Pi(0) = Pg(0) = 1- The Wilson coefficient C^ R for the b 
decay at m& is 



Cf R f (m b ) = 0.21 Ci R {m w ) - 0.17 \C u m (m w ) + C c 4R (m w )} 
-[0.37 PL(x t )- 0.015 P'{x t )]Cl R {m w ) . 



S7 



(42) 
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Let us now consider the experimental limits for the various coefficients in Eq. (f42|) . First 
of all, the R-parity violating couplings with n = 1 are strongly constrained by the non- 
observation of nucleon-antinucleon oscillation and double nucleon decay; A" 13 < 5 x 10~ 3 
and X" 12 < 10~ 6 with sfermion mass of 300 GeV |3l| . Thus C% R are negeligibly small. The 
constraint for Cf R {mw) comes again from the B — > K°7T° decay. We estimate the matrix 
element of the operator Of R as 

< K°n°\Oi R \W >^ i {m 2 B - ml) f K F^{m 2 K ) . (43) 

Using the similar values used in Eq. (^), we obtain 

|CJb| < 0.15 . (44) 

Note that this bound is also consistent with the data for the decay mode B + — > K°tt + [32 



Concerning the coefficient C% R , the consideration of the B — > J/tpKs decay gives |T4 

|C2h|< 0.02. (45) 
Applying again the bound ( f44|) to each component of Cf R , we get the bound on the C\ R ; 



2 

771 ~ 

\Cln\ < 015 -f. (46) 

Thus, to get \Cj R \ > 0.2 leading to a nearly maximal mixing-induced CP asymmetry, we 
need a sizable C\ R which can come about for > 3.5 mg c with x t — 1 

We conclude this section by pointing out that the direct CP asymmetry in the decay 
B^ — > tx^K can arise maximally as is the case with the coupling A^ 22 A^ 3 2 JTSfl. 



V. CONCLUSION 

We have discussed the effects of R-parity violation on the CP asymmetries in radiative 
B decays. When we allow R-parity and lepton number violating couplings which generate 
at one-loop level the tau neutrino mass of order 10 keV, they can induce rather large CP 
violating effects in the b — *> S7 decay. The direct CP asymmetry can be as large as 17 % 
if the R-parity conserving supersymmetric contribution is comparable to the SM one. The 
mixing- induced CP asymmetry can be almost maximal depending on the sfermion masses. 
For these sizable CP violating effects, it is required to have moderate sfermion mass splittings 
by factor 4 or bigger. If the atmospheric neutrino data from the Super-Kamiokande are to 
be explained by the oscillation between the muon and tau neutrinos whose masses are very 
light m v }l , T ^ 1 eV, then the effects of R-parity violation on the radiative B decays are 



3 The constraint on the relevant single baryon number violating coupling comes from T(Z — > 
ll)/T(Z -» hadrons) H, which gives | ^312 313 323 1 ^ 0-5 for rh — 100 GeV. This constraint is so 
weak that |Cg R4R | ~ 0(1) is easily allowed. 
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negligible. A large mixing-induced CP asymmetry is also possible with the R-parity and 
baryon number violating couplings with the similar order of sfermion mass splitting. 

These results could be contrasted with the R-parity violating effects on the CP asym- 
metries in the hadronic B decays such as — > J/ijjK, <pK or B^ — > tt^K, which could be 
significant without sfermion mass splitting and are rather insensitive to the neutrino mass 
restriction. 



Note added: While our work was being prepared, we encountered the paper 34 which 
considers the R-parity violating effect on the radiative B decay. It also deals with the RG 
running of the enlarged set of the Wilson coefficients which overlaps partly with our paper, 
and we find discrepancies in anomalous dimension matrix elements and the relation like 
(A10) and (All). 
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APPENDIX A: RG EQUATIONS WITH THE EXTENDED OPERATORS 



Here, we present the anomalous dimension matrix for the 28 operators including the 
standard set Oi ... 8 and the extended set Of 6 with q = u,d, s, c, b. We drop the indices L, R 
for the left-handed and right-handed set of operators as they have the identical anomalous 
dimension matrix. Omitting the usual 8x8 matrix for the standard set of operators, we 
have the following nonvanishing block- diagonal elements of the whole 28 x 28 matrix. 

The 2x2 submatrices mixing the operators Of 4 and 0\ 6 with themselves are 

-2 6 \ , (2 -6 \ /Ai . 

and ( \, (Al) 



6 -2 J V0 -16 

respectively. The 2x8 submatrix mixing the operators Of 4 with the standard set Oi,...,8 is 
given by 



fl fl _2e 2r _2r 2r _232_A ^ _l_ ZQA 

u u a u ab,s i u ab,s a u ab,s q w oii,s ai u ab,s J ~r ov o, 



a"qb,s r)Vqb,s n'Jqb,s •z'Jqb,s Vqb,s <-> T 27 "qb,s 

2 2 _2 2 416 r _ 232 r 70 , q r 

9 3 9 3 81 °1 U 81 °V 27 ou <lb,s 



(A2) 



where 



^ _ ( 1 for q = b or s 
qb)S \ otherwise, 

r <r [ 1 when q is an up-type, or down-type quark 

"qui Oqd 1 ~ 



otherwise. 

Similarly, the 2x8 submatrix mixing the operators 0\ 6 with the 8 standard operators is 
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32 A 
-gOqb 



■3 — Y^gfo 



2 2 _2 2 _44§A _|_ 200 f 119 u 

9 3 9 3 81 °1' a 81 0< ? d 27 



(A3) 



where 



1 for q = b 
otherwise. 



With this anomalous dimension matrix, the low energy Wilson coefficients are given by 

C e ff(pi) = U°"(n,n w )&"(n w ) (A4) 



where 

U eff (fi,fi w ) = V 
where V diagonalizes j e ff T 



iff / ' 



V 



-l 



(A5) 



D 



(A6) 



and is the vector containing the diagonal elements of the diagonal matrix of eigenvalues 
of 7 e// . 

By solving Eqs. (A4)-(A6), we find the following low energy Wilson coefficients at fib in 
terms of nonzero Cq^/jLw), Cq' c (/iw) in the L or R sector induced at fiw by R parity violating 
coupling X' n2j X'* 3j or X' nj2 X'* j3 , respectively, and C 3 f c {fi W ) in the R sector by A£ 12 A£ 3 , and 
the SM contribution C 2 l{^w) in L sector: 

o 10 

c 7 e/ V) = v 16/23 c e 7 ff (vw) + 1 (v u/23 - v 16/23 ) c e 8 ff (v w ) + c 2 { m ) E 



10 



10 



i=l 



10 



1=1 



+ (c^ w ) + ct{vw)) E r lv ai + c 6 W) E r-v ai + (c;W + QGM) E nV 4 

1=1 

10 10 

+ C^ w ) E rfv ai + C?(^) E r?i/* 

i=l i=l 
10 

+ + C 3 c (/i W ))E^V l + (QW) + Cl^w)) E^V 4 (A7) 

i=l 

10 

C 8 e// (/i fe ) = rj u ^C e 8 ff (f, w ) + C 2 ( m ) E M ai 

i=l 

10 

+ (c^ w ) + E nV* + c 6 W) e nV* + (c;w + E nV 4 

i=l i=l i=l 

10 10 

+ cf (nw) E nV* + ctivw) E nV* 

i=l i=l 

10 10 

+ (c 3 u (^) + c^ w )) e ffi/** + (c u 4 (fi w ) + ctQiw)) E nV* (as) 

i=i i=i 
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10 
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CM = ^(v 6/23 + v- 12/23 )c 2 M 
ci^\ii b ) = o 

cr^ b ) = I <y /23 + v~ i2/23 ) <#*w) + \ (v 6m - v- i2/23 ) cT d ' c M 

C^ b ) = I {v 6/23 + V- 12/23 ) Cl d ^ w ) + \ (^ 23 - r 12/23 ) C^w) (A9) 

where the quantities a i: hi, h i: rf , ff are shown in the TABLE I. 

The explicit numerical expressions with the choice nw — %, fJ>b — 4.8 GeV, a s (/iw) — 
0.120 and a s {n b ) = 0.214 are 

C% f (ji b ) = O.QQ87C 7L (fi w ) + 0.0920C 8L (fi W ) - 0.1732C 2L (/i W ) 

- 0.0974 (C* L (fi w ) + C s 6L (fi W )) - (0.6687 + 0.0875)C 6 fe L (M 

C 8 e {V) = 0.7032C 8L (M - 0.0801C 2L (M 

+ 0.1893 (C^w) + C s 6L (fi w )) + 0.3670C 6 6 L GM (A10) 
C$(lib) = 0.6687C 7il (/^) + 0.0920CW ? (/i w ) 

- 0.0974 (C* R (fi w ) + C s m (fi w )) - (0.6687 + 0.087 '5) C b 6R (»w) 

+ 0.2506 (Q R (fji w ) + C$ R {nw)) - 0M70C d R (fx w ) + 0.0880C^(^ W ) 
+ 0.0147 (C% R (»w) + C c m (fi w )) - 0.1732 (C 4 U R GM + C c m (^w)) 
CtH^b) = 0.7032C 8R {fiw) 

+ 0.1893 (C$ R {n w ) + C| fl (^)) + 0.3670^/^) 
+ 0.1893 (Q R (fi W ) + C c 6R (fi W )) 

- 0.0894 (C&Omw) + CU^w) + C c m (fi w )) 

- 0.0801 (C d 4R (fiw) + Ct R ^w) + C c 4R (fi W )) (All) 

where we have used ((J>w) — C 7 (fiw) — Cg(fiw) an d Cg (/ivy) — Cs^wO 011 the right- 
hand sides of the above equations. 
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TABLES 
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TABLE I. The magic numbers with R-parity violation 



X 


P 7 (x) 


Pg{*) 








0.197 


0.121 


0.157 


-0.0558 


1 


5/14 


1/4 


3/10 





/W 

1,200 J 


0.414 


0.301 


0.354 


0.0378 


(170\ 2 
\300 J 
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1 


1 


1 


1 



TABLE II. Values for the loop functions defined in the main text for several x's. 
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FIGURES 




FIG. 1. One-loop diagrams generating dipole moment interactions for the b — > 57 decay: the 
left-handed (a) and right-handed type (b) from the lepton number violating couplings, and the 
right-handed type (c) from the baryon number violating couplings. 
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FIG. 2. Scattered plots for Aqp as a function of B(B — > X s j) x 10 4 for the corresponding 
values of i]. The SM prediction is marked as a filled square. The numbers are values of \C\ L \ for 
the corresponding contour lines. The values of parameters of each case are shown in the text. 
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